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In this paper, we present a mathematical model for the angular projection of a rectangular 
arrangement of points in a grid. This simple, yet interesting problem, has both a scholarly value 
and applications for data extraction techniques to study the physics of various systems. Our work 
can interest undergraduate students to understand subtle points in the angular projection of a grid 
and describes various quantities of interest in the projection with completeness and sufficient rigour. 
We show that for certain angular ranges, the projection has non-distinctness, and calculate the 
details of such angles, and correspondingly, the number of distinct points and the total projected 
length. We focus on interesting trends obtained for the projected length of the grid elements and 
present a simple application of the model to determine the geometry of an unknown grid whose 
spatial extensions are known, using measurement of the grid projection at two angles only. Towards 
the end, our model is shown to have potential applications in various branches of physical sciences 
including crystallography, astrophysics and bulk properties of materials. 


I. INTRODUCTION 

Many physical systems have order in the arrangement 
of their constituent elements - the study of which enables 
us to understand the structure and various properties of 
the system. A sound understanding of this internal order 
can help us devise relevant data extraction techniques 
which harness prior information about the order to help 
probe the system better. A typical example of utilizing 
the order for the study of the system can be directly seen 
in, but is in no way limited to, crystallogr aphyUlU , such 
as Bravais lattices Q, where physical properties of the 
system such as the material density, electronic properties, 
Young’s modulus, Bulk modulus, coefficient of thermal 
expansion etc. are dependent on the crystal structure 
(order) of the lattice. 

In this paper, we study a simplistic arrangement of 
a rectangular grid of constituent elements; which is a 
common, two dimensional geometrical arrangement seen 
in many physical scenarios. The grid elements are ar¬ 
ranged in a rectangular fashion with constant spacing 
between neighbouring elements and this constant spac¬ 
ing or the period, may have different values along the 
two perpendicular directions of the rectangular grid. We 
mathematically model the angular projection of the grid 
elements onto a circular screen and try to understand the 
angular dependence of the projection pattern on various 
grid parameters. In particular, we study the variation of 
the projected length, the number and distribution of pro¬ 
jected points with the angle of projection, and try to cor¬ 
relate the modeled trends with the grid parameters. The 
constructed model will be shown to have practical im¬ 
plementations in physics; as data extraction techniques 
- which combine a physical model of the structure with 
observations to make inferences about the system proper¬ 
ties. As we will show, the formulated model can be used 
to determine the number of grid elements, the size of 


each element and spacing between them in an unknown 
grid, using only two angular projection measurements. 
With a simple setup coupled with tools and methodology 
available to young students, they can appreciate both the 
reflective and application-oriented nature of this work. 

In section im we describe the grid setup, describing 
various parameters in detail and explain the idea of the 
angular projection of the grid. In section UTTl we build the 
mathematical model of the angular projection and calcu¬ 
late quantities of physical interest which is extended in 
section HVl where we discuss a methodology based on the 
model to determine the grid parameters using a two-angle 
measurement approach. In Section [Vl we outline a num¬ 
ber of possible applications of our work and its possible 
extensions in various branches of physical sciences such as 
crystallography, astrophysics, signal encoding and bulk 
properties of mechanical systems. We conclude in sec¬ 
tion IVII by summarizing our work. 


II. THE GRID SETUP 

Consider a rectangular grid of points having n rows 
and to columns, that is a total number of No = nm 
grid points. Each grid point has an associated circular 
grid element of radius ao- For convenience, we have taken 
both n and m to be odd positive integers and have chosen 
an orthogonal Cartesian coordinate system having axes x 
and y, and the origin is coincident with the center point of 
the grid. Immediately, we see that the grid is symmetric 
with respect to both x and y- axes. Let Xq and yg be the 
constant spacing between neighbouring grid point centers 
along the x and y directions, respectively. The grid setup 
is shown in Figure ©• We use discrete Latin indices (i, j) 
to label the grid point center in the i-th column from the 
right and j-th row from the bottom (ref. Figure© such 
that a general point P on the grid can be labeled as, 

P = P i:j = P{xi,yj) = (y(m —2* + 1), y(n- 2j + 1)) , 

(1) 
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FIG. 1. An Example Setup of the Rectangular Grid of 7 x 
5 points, each associated with a finite sized circular element 
(size exaggerated for clarity i.e. not to scale). Also shown 
are the spacing between neighbouring grid centers, a general 
projection angle, total projected chord length and the circular 
screen. 

with i G {1, 2, 3,..., to} and j G {1, 2, 3,..., n}. The grid 
is, therefore, completely characterized by n , to , Xq, yo and 
ao, which we call the parameters of the grid. The total 
grid length along the x and y-directions is, respectively, 
L x and L y : 

L x = (to - l)a; 0 ; L v = (n - l)y 0 (2) 

The angular projection of the grid at an angle 9 (as shown 
in Figure (JTJ ) is the projection of all TVo grid elements 
onto a circular screen of radius R, also having its cen¬ 
ter at the origin of the coordinate system. The angle 
9 is measured from the negative y-axis with 9 G [0, 27r] 
and is conventionally taken to be positive in the counter¬ 
clockwise direction. We exclusively work in the range 
9 G [0, f-], since for higher angles, the extension of the 
results is straightforward due to the two-fold symmetry 
in the grid setup. This projection can be physically real¬ 
ized by passing a light beam (with a plane wave-front )[|| 
through the grid at different angles, and detecting it on 
the concentric circular screen screen kept in the x — y 
plane. The shadow regions, against an illuminated back¬ 
ground on the screen will be the projection of the grid el¬ 
ements. These projections can be detected using a photo¬ 
detector Q which can move along the screen to measure 
the changes in intensity of light due to the grid pattern. 
With the grid setup so ready, we now model the angular 
projection and relate it to the grid parameters. 


III. MODELING THE ANGULAR PROJECTION 

We now model the angular projections of the grid onto 
the concentric circular screen and develop ways to utilize 


FIG. 2. Coverage Fraction (solid line) and Total Projected 
Chord length (dashed line) of the Angular Projection of a 
Grid having n = 5, m = 7 and xo = 1, yo = 2, R = 10 and 
ao = 0.07 in standard distance units. 

this model in physical situations for appropriate data ex¬ 
traction applications. We first consider the extent of the 
projection on the circular screen, that is, the length be¬ 
tween the two most separated (projected) points. It is 
readily seen that this length 1(9) will be the projected 
length of the grid diagonal. We first get an expression 
for the projected length of the intercepting chord, which 
can then be used to find the arc length of the projection 
along the screen. The chord length can be understood as 
the projection of the line joining the diagonally opposite 
grid centers and an additional projected length of 2ao 
due to the finite radii of the grid elements located at the 
diagonally opposite ends, 

1(0) = L x cos 9 + L y sin 9 + 2ao . (3) 

This can be seen to have a maximum value for 9 = 
tan ~ 1 (L y /L x ). The length of the projected arc along the 

circular screen will be denoted be 1(9) and can be found 
using the fact that the angle subtended by the projection 
at the grid center is [2 sin - 1 (l(9)/2R)] , where R is the 
radius of the screen, such that, 

W) = 2-Rsin -1 • (4) 

Naturally, in the angular range [0, 27 t] the projected arc 
length of the grid is an oscillatory function of the projec¬ 
tion angle 9. In our range [0,7r/2], the projected chord 
length of the grid is shown as the dotted line in Fig¬ 
ure [2] which corresponds to a grid having n = 5, to = 7 
with Xo = 1, yo = 2, R = 10 and ao = 0.07 in stan¬ 
dard distance units. We now go on to determine the 
angular dependence of the number of distinct projected 
grid points, the projected length due to the grid elements 
and other related parameters useful for a data extraction 
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model. At any general 9, we expect a maximum num¬ 
ber of N 0 projections of the grid elements, each having a 
one-dimensional projected length of 2ao on the projection 
chord. In certain angular ranges, however, the projection 
of some grid elements may overlap with one another, and 
in such a case, each grid element will not project itself 
distinctly. We call this the degenerate scenario, since 
multiple elements on the grid will have a finite overlap 
with other grid elements. At certain specific values of 0, 
there would be a complete overlap between different grid 
elements and at such angles, the number of projected grid 
elements will decrease. Around each such specific angle, 
there would be a finite angular width, within which even 
though the grid centers will be distinct, there will be 
overlap between different grid elements due to the finite 
extent (~ ao) of each element. But the projection at 
a non-degenerate angle will have all the Nq = nm grid 
elements (and therefore, the grid centers) projected dis¬ 
tinctly, each spanning 2ao in length. Let Lp(9) be the 
projected length of the chord due to the grid elements at 
angle 6. For non-degenerate angles, the projection due 
to the grid elements will have a length given by, 

Lp(9) = 2nmao , 9 is a non-degenerate angle . (5) 

However, this condition of non-degeneracy must be ac¬ 
companied by a constraint on the size of each grid ele¬ 
ment. We need to ensure that at non-degenerate angles, 
there is no overlap between grid elements due to a large 
ao since all element projections are confined to a finite 
chord length 1(9). We impose the following strict condi¬ 
tion on the size of each grid element in view of this, 

( Lp(9) = 2nmao) < min(Z(0)) . (6) 


This condition essentially demands that the projected 
length of the grid elements L P (9 ) should not exceed the 
total expanse of the grid projection 1(9). We make it 
strict by imposing that the maximum possible length of 
the grid element projections must be smaller than the 
smallest value of 1(9), ensuring there would be no over¬ 
lap at any non-degenerate angle due to the finite size of 
each element. Though it must be pointed out that the 
minimum of 1(9) could be either L x + 2ao or L y + 2ao, 
and correspondingly Lp = 2mao or 2?iao, respectively, 
and not 2nmao. Thus, the condition of Eq. © presents 
a very strict limit on ao which can be succinctly written 
as, 


ao < 


L X + Ly 

A(nm — 1) 


(7) 


Please note that the grid element size a o is exaggerated 
in Figure [T| for visual clarity. The grid element radius 
ao must follow Eq. 0. The next crucial aspect is to 
understand the degeneracy in the grid projection due to 
the angular orientation of the projection i.e. when their 
projections start overlapping at certain angles of projec¬ 
tion. Let N(9) be the number of projections of distinct 
grid element centers at angle 9. It is easily seen that for 


all those angles such that 9 ^ tan 1 (pxo/qyo ), with any 
combination of p , q satisfying, 

P G {0,1,2,..., (m — 1)} , (8) 

q G {0, 1,2, ...,(n-l)} , 

p = q = 0 is not allowed , 

the number of distinct centers is N(9) = No = nm, since 
for these angles, all the grid centers will have a distinct 
projection. We label such angles by the set {a}, such 
that whenever 9 G {a}, each grid center has a distinct 
projection. Let, x = t&n~ 1 (pxo/qyo) be a subset of all 
possible angles 9, labeled by the set {x}, for which all 
grid element centers do not have a distinct projection. 
Thus, whenever 9 G {\}, there exist degeneracies in the 
projection of the grid and this is termed as the “full- 
eclipse” scenario, since the alignment between the grid 
elements is perfect. In the full-eclipse case, we have, 
as expected, N(x) < N 0 . The number of such angles 
in the set {%} is the number of unique combinations of 
(p/q) in the lowest rational form satisfying Eq. ©. As 
an illustration of this point, we see that p = 0 for any 
q / 0 corresponds to x = 9 = 0, in which case, there 
are only m distinct projections of the grid centers, one 
due to each column. On the other extreme, q = 0 for 
any p / 0, will correspond tox = 0 = 7r/2, and there 
would be only n distinct projected centers, one due to 
each row. We label each angle in the set {x} with a 
subscript “k” such that a full-eclipse angle is labeled as 
Xk- At Xk, many grid elements will be degenerate be¬ 
hind a single grid element and there will be a total of 
N( Xk ) number of distinct grid elements projected on the 
screen. Around each Xk > there will be an associated an¬ 
gular width A Xk associated with it, within which there 
will be overlap between projection of grid elements, even 
though the element centers are distinct. In the range 
{Xk ± Ay,,.} i.e. when 9 G [xfc - A Xk, Xk + A Xfe], the el¬ 
ement centers will be distinct but the projection of the 
grid elements will overlap such that each of the N(xk) 
elements which were distinct at Xfc will now be “smeared 
out” due to the finite size of each element. We first cal¬ 
culate the number of distinct grid centers projected at 
a given degenerate angle Xfc and the projected length of 
the grid elements, before moving on to the more involved 
case when 9 G {xfc ± Ay/c}. Only rational combinations 
of (p/q) in the lowest form satisfying Eq. © will con¬ 
tribute to a unique X- To ensure that we work only with 
unique angles x while considering all combinations of Eq. 
©, we reduce (p/q) to its lowest rational form and write 
it as (u/v), such that u and v are co-prime. Thus, there 
is a many-to-one map between the pair (p, q) to the (u, v) 
pair, such that each (p, q) combination can be mapped 
to a corresponding ( u , v) pair, and hence, to the angle 
X = t&n~ 1 (uxo/vyo). It is worth mentioning here that 
the lowest rational combination for p = 0 for any q ^ 0 
corresponds to the (0,1) ( u , v pair and for q = 0 for any 
p / 0, to the pair (1,0). Let, for a given (u,v) pair, 
the total number of points (grid centers) which have a 
non-distinct projection be j3(u, v) = /3(xk), such that the 
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number of distinct projected points will be, 

N(Xk) = N 0 -P(xk) = nm-Pixk) ■ ( 9 ) 

As mentioned in section HT1 we label the row of grid points 
closest to the bottom of the grid as “Row-1”, the next as 
“Row-2” and so on till “ Row-n ”. We analyze the rows for 
the projection of their points in this numerical sequence 
to determine the total number of distinct projections. 
The number of points of Row-t having a non-distinct pro¬ 
jection due to overlap with projections of other points in 
the same row or points of previously analyzed rows is 
given by j3 t (u, v) such that the total number of such de¬ 
generate (non-distinct) points in the grid is, 

n 

p(u,v) = Y^0t(u,v) . (10) 

t=l 

It can be seen with some thought that for Row-1, there 
are (m — u) points having a non-distinct projection if 
and only if v = 0 (which corresponds to u = 1 and 
X = §), else all points will have a distinct projection. 
Thus, when (u,v) = (1,0), we see that (in — 1) points 
of Row-1 have a degenerate projection, and there is only 
u = 1 distinct projected point for this row. From Row-2, 
there would be u points having distinct projections, so 
that (m — u ) of its m points have a non-distinct projec¬ 
tion only if v = 1 as they would overlap with the last 
(to — u) points of Row-1. We emphasize that the readers 
give it sufficient thought to convince themselves that for 
Row-t, there would again be u distinct projections and 
correspondingly (to — u) of its points having non-distinct 
projections only if v G {0,1, 2,3,..., t — 1}. The value of 
v determines with which row the (m — u ) points of Row-t 
will overlap, such that v = 0 corresponds to the row it¬ 
self, v = 1 signifies overlap with the points of Row-(t — 1) 
and so on. This can be expressed by formulating an ex¬ 
pression for f5 t (u,v ) with the use of Kronecker delta]!} 
functions as follows, 

t—i 

Pt(u, v) = (to - u) S v> k • (11) 

k =o 

With an expression for (3t(u, v) ready, we can easily find 
the total number of grid points which have a non-distinct 
projection using Eq. m, 

n t—1 

p(u,v) = - u)S v , b , (12) 

6=1 6=0 

which can be further simplified as, 

n— 1 

P(u, v) = (to - u) y~^(n - s)6 Vt „. (13) 

s=0 

For any given y, v £ {0,1, 2,..., (n — 1)} from Eq. 
depending on the values of p, q and u. Thus, the Kro¬ 
necker delta 6 VtS selects only one term in the sum with 
v = s and we get an expression for /3(u,v) as, 

/ 3(u,v ) = (to — u)(n — v) . (14) 


The total number of distinct projected grid element cen¬ 
ters for a given Xk can be calculated from Eq. ©, 

N ( u, v ) = N (xk) — nm — (m — u)(n — v) = mv + nu — uv. 

(15) 

This can be made more suggestive by writing in a form 
which allows comparison with Nq , 

N(xk) = N 0 (- + --^-) . (16) 

\n to A o / 

Each of these N(xk) distinct points may have a finite 
number of other grid elements exactly eclipsed behind 
these distinct elements. The special cases of 9 = 0 and 
9 = ^ can be easily recovered for which N is m and n, re¬ 
spectively. We can see that N(xk ) < -W) since the factor 
in the parentheses in Eq. (USD is less than unity. This is 
expected because some elements will have an overlapping 
projection and the total number of distinct projections 
will decrease, as is the essence of the full-eclipse case. For 
these degenerate angles Xk, the total length of projection 
due to the grid elements will be denoted by Lp(xk) and 
given by, 

Lp( Xk ) = 2N( Xk )a 0 , (17) 

since each of the N(xk ) points contributes 2ao to the pro¬ 
jection length of the chord. As we change 9 in the range 
{Xk ± A^t}, the full-eclipse degeneracy breaks and the 
grid elements behind every distinct element will begin to 
separate until all the No = nm elements have completely 
distinct projections. In this angular range, there will be a 
finite overlap between the projections of the grid elements 
which were degenerate at Xk and this would happen for 
all such N(xk) points. This angular width and the corre¬ 
sponding projection length due to the grid elements will, 
in general, depend on the number of overlapping elements 
with the concerned distinct element at at Xk- Using Eq. 
m, we see that “Row - t n has, 

t -1 

Qt(Xfc) = m - Pt(u, y) = TO - (to - u) S v> k , (18) 

k=0 

number of distinct element centers, which can be labeled 
as {1, 2,3,.., fit(Xk)} from the right. Thus, each of the 
N(xk ) distinct element centers at 9 = Xk can be iden¬ 
tified with a coordinate ( t,u(t )) with 1 < t < n and 
1 < oj(t) < f \(xk) such that t labels the row and uj(t) 
the column. In the full-eclipse scenario, it can readily 
be seen that for the distinct grid element (t,u>(t)) will 
be overlapping with a total of J(t,oj(t),Xk) number of 
elements at Xfc, given by, 

T . . . . . f,m — uj(t), ,n — t,\ _ 

J(t,u(t),Xk) = min -— J, L -J + 1 , (19) 

\ U V J 

where [-J is the floor function^ (The additional unity 
comes from the element itself). We now evaluate an ex¬ 
pression for A Xk-, the angular width around Xk needed 
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to be swept to project all J(t,ui(t),Xk) points separately. 
It can be seen that this width will only depend on (u, v) 
and ao, and not on J(t,uj(t),Xk) since the projection is 
done using a plane wave-front. We define the distance 
between the centers of two adjacent, overlapping points 
at Xk as n which has the form, 

9 = ^u 2 x% + v 2 yl . (20) 

At 9 = Xk ± Ayfc, the projection of all the J(t,u)(t),Xk ) 
overlapping points, for each of the N(xk) points will be 
just separated and tangential to each other, such that, 

A Xk = sin -1 ■ ( 21 ) 

At an intermediate angle 9 G [xk ~ &Xk, Xk + Ax*], 
there will be some finite overlap between grid elements 
that were degenerate at Xk, for each of the N(xk) points. 
The projected chord length of the grid elements due to 
the overlapping points associated with the (t,ui(t)) point 
can be written as, 

Lp,t(0, Xk) = [J(t,uj(t),Xk) - l]^sin|0-Xfc|+2a o , (22) 

where the subscript “t” denotes that the point (t,w(t)) 
is under consideration and the argument Xk refers that 
the angle is being varied around 9 = Xk in the range 
[Xk - A Xk,Xk + A Xk]- For all N(xk) such sets of over¬ 
lapping points, the total projected grid element length 
can be obtained by summing over both t and u>(t) in Eq. 

m , 


Lp(9,\k) 


n Cl(t) 

E E 1) 


t= 1 w(t) = 1 


fism\9 -Xk\ 


+ (2ao) 


i n(t) 

E Ed) 

£—1 oj(t) = 1 


(23) 


It can be identified that the second double summation in 
Eq. (E3l) is the total number of distinct grid elements at 
9 = Xk and therefore, Eq. (E^l) may be written as, 


Lp{9,Xk ) 


n 0 (£) 

E E (■£(*.^(t),xk) -1) 


£—1 u)(t )=1 


+ 2aoN(xk) ■ 


Msin 1 9 - Xk I 
(24) 


We now define a quantity called the “coverage fraction” 
which is a dimensionless number, giving us the fraction 
of the total projected chord length of the grid 1(9) which 
is covered by the projection of the grid elements. On 
the circular screen, this will correspond to the fraction 
of the length over which the photo-detector will detect 
a lower intensity than the background illumination due 
to the source plane wave-front which passes through the 
grid. This quantity A (9) can be measured and inferences 


about the grid can be made from these observations, as 
will be elaborated in the next section. It can be defined 
as follows, 


m = 


Lp(9) 

1(9) 


(25) 


where 1(9) can be used from Eq. (0 and Lp(9) 
from Eq. (l24l) or ©, depending on whether 9 G 
[Xfc — Axfc,Xfe + Axfc] or not, respectively. Whenever 
6 G [xfc ^ Axfc, Xk + Axfc], the coverage fraction would 
be lower, owing to the overlap in the grid projection. 
Thus, for each 9 = Xk, there would be a finite angular 
width in which the coverage fraction would decrease from 
its non-degenerate value. In the next section, we present 
an application of our model of the grid projection to de¬ 
termine the unknown grid parameters n, to, xq, yo and ao 
of a grid whose spatial extent is known. 


IV. GRID PARAMETER DETERMINATION - 
THE TWO ANGLE APPROACH 


With a mathematical model of the grid projection 
ready, we present a formulation to determine unknown 
parameters of the grid whose spatial extent is known. 
This will be done by measuring of the first two angles 
at which degeneracy exists in the projection. Consider a 
situation in which we have a grid with a known extent 
in x and y directions i.e. L x and L y are known. The 
grid parameters n,m,Xo,yo and a o are unknown and we 
wish to find these, by measurement of the two smallest 
non-zero angles xi and X 2 at which there is degeneracy 
in the grid projection. Physically, this can be achieved 
by continually changing the angle of the light beam (as 
described in section ED from 9 = 0 to increasing values 
and measuring the projected length of the grid elements 
on the screen, xi is the first or smallest full-eclipse angle 
at which there is a fall of the grid element length pro¬ 
jection and X 2 will be the next largest such angle. With 
the two angles known, we would be able to determine the 
grid parameters as demonstrated here. It is easily seen 
that these two angles will be such that they follow, 

tanyi = j— , (26) 

L/y 

tanx 2 = —^—■ (27) 

Ty 2/0 

Once we have measured Xi ; we can readily infer xo using 
Eq. (1261) since L y is known. We can then use Eq. CD and 
our knowledge of L x to determine the number of columns 
in the grid, 


I.t 


L y tan xi 


1 . 


(28) 


The measurement of the other angle X 2 can be then used 
to deduce yo using Eq. m to get, 


Vo 


xq 

V tan X 2 ’ 


(29) 
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which can be then used together with Eq. dm to get the 
value of the number of rows n in the grid, 


tan x 2 

n = - 

tan X 2 - tan xi 


+ 1 . 


(30) 


To determine the grid element size ao, we would need the 
value of angular width where the first degeneracy occurs 
i.e. 2Axi which can be measured by the photo-detector. 

At xi) we see that fi = fix = + (n — l) 2 £/o) and 

N (xi) = m(n — 1) + 1. Substituting in the expression for 
Axi at 9 = xi from Eq. (l^Tl) . we get an expression for 
ao as, 


ao = -Mi sin(2Axi) . 


(31) 


Thus, with the knowledge of the physical extension of 
the grid and measurement of the two smallest angles for 
the degenerate case, we are successfully able to deduce 
the grid parameters using a suitable application of the 
mathematical model developed in the previous sections. 


V. DATA EXTRACTION APPLICATIONS TO 
MODERN SCIENCE 

In this section, we present a short discussion on possi¬ 
ble applications of our model to different setups in mod¬ 
ern science. The model, which uses the projections of a 
rectangular grid onto a circular screen, in essence pro¬ 
vides us with a methodology for data extraction tech¬ 
niques. The occurrence of an ordered structure in physics 
and other branches of science is fairly common and in¬ 
dividual problems can motivate specific adaptations of 
the model. A given grid structure in three dimensions 
can be thought of consisting of many two-dimensional 
slices of grids and our model can be applied to each of 
these slices separately, while also be used to estimate the 
spacing between adjacent slices. (Ref. section HVl) . The 
model harnesses the symmetry in the structure, which in 
our case of a rectangular grid is a C 2 symmetry about the 
x and y axes. The total projection of the grid onto the 
screen 1(9) will be representative of this symmetry if 9 is 
allowed to vary from 0 to 27r. A simple plot of 1(9) with 
9 will easily give us estimates of the spatial extent of the 
structure i.e. L x and L y in our case. For grid arrange¬ 
ments having a different symmetry such as triangular, 
hexagonal or circular etc.; the symmetry properties[1, Hj 

can be obtained by studying the polar plot of 1(9) for 
9 e [0,27 t]. The application of the model to classical 
macroscopic systems can be thought which use the an¬ 
gular pattern of the grid projection to infer properties 
of the system. Characterization of periodically arranged 
mechanical systems (eg. 0) can be implemented using 
this model. Our model can also be used to study the 
physical properties of a material, in particular its elastic¬ 
ity and response to temperature change. A known grid 
whose angular projections are known can be placed on 


a given material, which is then subject to changes. An 
external stress would tend to produce a strain in the ma¬ 
terial, altering the spacing between the grid points and 
may even alter the size of each grid point. A similar ef¬ 
fect would be produced by the application of an external 
pressure. By measuring the changes in the grid projec¬ 
tion pattern due to these external influences, one can 
determine elasticity [ll] properties of the material, such 
as its Young’s Modulus flD] and Bulk Modulus ll^. This 
method can be easily extended to obtain the temperature 
expansion coefficient [l2j of the material since we expect 
changes in its physical dimensions due to a temperature 
change. Exact implementation of these techniques using 
the grid projection would need a detailed construction 
of the setup and a proper analysis of all underlying as¬ 
sumptions, if any - such as uniformity of properties, ac¬ 
curacy of the projection etc. We also motivate the use 
of this model in astrophysics, where it can be used to 
measure the curvature of space-time and for the detec¬ 
tion of gravitational waves. The model relies on changes 
in the physical extent of the grid due to variations in 
the local space-time, a characteristic common to ongo¬ 
ing large interferometric projects in this field like LI GO 
and LISA\ 13j. Application to signal encoding can also be 
developed, whereby the grid pattern can be the encoded 
signal which can be transmitted at the speed of light to a 
receiver. The key to the decryption are the grid parame¬ 
ters which are only known to the sender and receiver. The 
information is contained in the particular angle at which 
the projection is transmitted. Such a methodology can 
be used for high-speed secure data transfer. One of the 
areas of modern physics which deals with system having 
a regular structure is crystallography. Our method can 
serve as an addition to standard X-Ray Diffractionl 13- 
0 and Neutron Diffraction jl2| techniques to determine 
the underlying crystal structure. Using our model, we 
can infer the orientation of the lattice planes and also 
estimate the spacing between adjacent parallel planes. 
Proper extensions of this work can also be used to char¬ 
acterize a Bravais Lattice completely by measuring all 
six of its lattice constants. Further thoughts on appli¬ 
cations in crystallography can deal with measurement of 
the material density and detection and analysis of lattice 
defects 0 ]. We emphasize that angular projection is an 
important and useful methodology to study the proper¬ 
ties of a system since it is easy to implement and helps 
reducing the degrees of freedom of the data. In par¬ 
ticular, the projection of a two-dimensional grid can be 
taken onto a one-dimensional screen, thereby reducing 
a degree of freedom to work with. Our work should be 
seen as a framework for a data extraction model which 
can be adapted suitably for different physical systems to 
appropriately extract details of the underlying structure 
by measuring its angular projection. In addition to data 
extraction uses, undergraduate students can pursue this 
as an engaging exercise to understand interesting aspects 
of the projection of a regular grid. 
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VI. CONCLUSION 

We have presented a simple mathematical model for 
the projection of a rectangular grid onto a circular screen 
and discussed various applications to physical systems as 
data extraction techniques. Regular grids are a common 
geometric arrangement and our work tries to model the 
angular projection of the grid and relates the observed 
trend of physical observable quantities like the projected 
length, the number of projected points etc. with the grid 
parameters. We take into account the non-distinctness 
in the projection due to the overlap of projected grid 
elements in certain angular ranges, and give a sound 
methodology to deal with such cases leading to interest¬ 
ing trends for the coverage fraction. We presented an ap¬ 
plication of the model to determine the parameters of an 
unknown grid using the measurement of the first two an¬ 
gles at which degeneracy occurs. Our model can be seen 
to have relevance and uses to formulate data extraction 
techniques. Applications of the model are motivated in 
different branches of modern science such as macroscopic 


elasticity property determination, astrophysics, crystal¬ 
lography and signal encoding. Apart from such uses, 
such a model has a scholarary value to it. It can be 
of particular interest to undergraduate students who like 
working on simple, yet interesting, innovative problems 
to understand the structure of the solutions with suffi¬ 
cient rigour. 
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